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ABSTRACT 

The Naval Postgraduate School open circuit oscillating 
flow wind tunnel was used to study the blowing requirements 
to maintain an attached turbulent boundary layer in an oscil- 
lating freestream flow with an adverse pressure gradient. 
Boundary layer separation was visualized through the use of 
tufts. Freestream flow oscillation frequency was found to 
have an effect on the flowing required to maintain attached 
flow. This frequency dependence exhibits characteristics 
which suggest resonant behavior. 



4 



TABLE OF CONTENTS 



I. INTRODUCTION 11 

II. BACKGROUND 12 

A. NON-STEADY FLOWS ^ 12 

1. Analytic Studies 12 

2. Experimental Studies 13 

B. BLOWING TO AVOID SEPARATION 14 

1. Experimental Studies 14 

2 . Analytic Studies 12 

C. APPLICATIONS 18 

III. DESCRIPTION OF APPARATUS 19 

A. OSCILLATING FLOW WIND TUNNEL 19 

1. General Description 19 

2. Rotating Shutter Valve 20 

3. Test Section 21 

B. MODEL 22 

C. INSTRUMENTATION 23 

1. Freestream Sensors 23 

2. Blowing Coefficient 24 

3. Visualization of Separation 25 

IV. TEST PROGRAM 38 

A. EXPERIMENTAL PROCEDURE 38 

B. DATA REDUCTION 39 

V. RESULTS 

VI. CONCLUSIONS 39 

VIII. RECOMMENDATIONS 30 



5 



COMPUTER PROGRAM 



61 



LIST OF REFERENCES 64 

INITIAL DISTRIBUTION LIST 67 



6 



i 



LIST OF FIGURES 



FIGURE 

111.1. 

111. 2. 
III. 3. 
Ill .4. 
III. 5. 
III. 6. 
III. 7. 

111. 8. 

111. 9. 

111. 10. 

111. 11. 

V. 1. 
V.2. 

V. 3. 
V.4. 

V.5 

V.6. 

V.7. 

V.8. 

V.9. 



Plan View of Wind Tunnel 

Photograph of Rotating Shutter Valve 

Photograph of Shutter Valve Drive Mechanism 

Photograph of Wind Tunnel Test Section Closed-- 

Photograph of Test Section Open with Model 

Photograph of Test Section and Control Console- 

Drawing of Model and Plate 

Drawing of Plenum and Blowing Slot 
Configuration 

Schematic of Test Setup 

Pressure Drop Through Orifice vs. Volume Flow 
of Air 

Photograph of Diverging Planes from Down 
Stream 

Photograph of Tufts Showing Corner and Side- 
wall Separation 

Blowing Coefficient Required vs. Frequency, 

R^ = 5.93x10® 

Blowing Coefficient Required vs. Frequency, 

Ry = 7.31x10^ 

Blowing Coefficient Required vs. Frequency, 

Ry = 8.53x10® 

Pressure Coefficient vs. Position 

Pressure Gradient vs. Frequency 

Velocity Perturbation vs. Frequency 

Ry = 5.9 3x10 

Velocity PeiturbaLion vs. Frequency 

Ry = 7.31x10® 

Perturbation vs. Frequency, Ry = 8.53x10® 



PAGE 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

50 

51 

52 

53 

54 

55 

56 

57 

58 



7 



LIST OF SYMBOLS 



D 



m 



X 



R 



h 

a 

m . 



I.. 

R,- 



Width of the blowing slot, feet 

Drag coefficient 

Lift coefficient 

Moment coefficient 

Pressure coefficient , (P-P^)/q 

Cp derivative with respect to x position along plate 

Volumetric coefficient of blown air, as defined in 
Equation II-l 

Momentum coefficient of blown air, as defined in 
Equation II-2 

Blowing coefficient required to maintain attached flov7 

Height of the blowing slot, feet 

Slant length of the diffuser section, feet 

Mass flow of blown air, slugs per second 

Pressure at exit section of diffuser, pounds per square 
foot 

Total pressure inside blown air supply duct, pounds 
per square foot 

Freestream static pressure, pounds per square foot 
Volume flow of blown air, cubic feet per second 
Freestream dynamic pressure, pounds per square foot 
Slot height Reynolds number, V^h/v 
Reynolds number, U^Jl/v 

Characteristic area of interest, area of blown wall 
of diffuser section for this study, square feet 

Total temperature inside blown air supply duct , degrees 
Rankine 



Freestream temperiaturc , degrees Rankine 



8 



Velocity at exit section of diffuser, feet per second 
Freestream mean velocity, feet per second 

V. Jet velocity of blown air assuming isentropic expansion 

^ from P and T to P and T, feet per second 
o o 

V.j^ Minimum jet velocity that will just prevent separation, 
^ feet per second 

X Slant distance along diffuser wall measured from inlet 

section, feet 

APj. (Pj,-P^) , the pressure rise through the diffuser, pounds 
per square foot 

AU^ Amplitude of perturbation in velocity at inlet to the 
diffuser, feet per second 

6 Displacement thickness of boundary layer at exit sec- 

tion, feet 

p Absolute viscosity, slugs per foot-second 

V Kinematic viscosity, feet squared per second 

p DensiLy, slugs per-' cubic foot 

w Frequency of oscillation, radians per second 



9 



ACKNOWLEDGEMENT 



The author wishes to express his indebtedness to Dr. 
Howard L. Power, Assistant Professor, Department of Aero- 
nautics, for his encouragement and thoughtful guidance 
during the present work. 

Grateful acknowledgement is also due Dr. James A. 

Miller for his initial helpful suggestions, and to Dr. T. 

H. Gawain for his help in developing parts of the analysis. 

The author also wishes to express special appreciation 
to the technical staff of the Department of Aeronautics, 
and in particular Messrs. Theodore B. Dunton and Robert A. 
Besel, for their cooperation and assistance in this 
investigation . 



10 



I. INTRODUCTION 



The study of viscous unsteady flows have become increas- 
ingly important in recent years. Typical examples of un- 
steady flows are found in helicopter blade aerodynamics and 
in gas turbine engine internal flows. The development of 
advanced design concepts in helicopters, ships, and propul- 
sion systems requires a basic understanding of the nature of 
unsteady flows . The complexity of the problem does not lend 
itself readily to either analytic or experimental study, 
even when limited to flow with laminar boundary layers. 
Furthermore, most of the problems of interest are characterized 
by turbulence and turbulent boundat'’v layers . which greatly 
increases the difficulty of the study. Regardless of the 
difficulties involved the field is of such importance as to 
V7arrant the expenditure of time and effort to widen the un- 
derstanding of its characteristics. 

The present study was to experimentally examine the ef- 
fect of different levels of blowing on turbulent boundary 
layer separation in an oscillating freestream. The effects 
of changes in frequency and magnitude of the freestream 
oscillation on the blowing requirements to maintain an at- 
tached boundary layer v;ere studied. 
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II. BACKGROUND 



A. NON-STEADY FLOWS 

1 . Analytic Studies 

Previous analytic studies have been confined mostly 
to flows which have a laminar boundary layer. The earliest 
was Stokes [Ref. 1] study of the doubly infinite flat plate 
oscillating in its own plane in a fluid at rest. Rayleigh 
[Ref. 2] considered the second order effects of this prob- 
lem. Schlichting [Ref. 3] expanded Rayleigh's work to in- 
clude the boundary layer assumptions. All of these were 
simple unsteady flows without a mean flow or pressure gra- 
dient. Lighthill [Ref. 4] later treated the case of sm.all 
magnitude, low frequency, sinusoidal oscillations super- 
imposed on a steady flow about a cylinder and a semi-infinite 
flat plate. Using small perturbation theory and resti'>icting 
the solution to first order terms he determined that, at low 
frequencies the flow is essentially quasi-steady. The un- 
steady flow had the same characteristics for any instantan- 
eous magnitude of freestream velocity as a steady flow of 
the same velocity. As frequency is increased a limiting 
frequency is reached beyond which the boundary layer reacts 
as it would without a mean flow. Lin [Ref. 5] analyzed the 
high frequency case for large amplitude oscillations. He 
found that the governing equations became essentially lin- 
ear at high frequencies allowing separation of oscillatory 
and steady components. Nicker'son [Ref. 6] expanded this 
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work to include higher order perturbation terms. Barriol 
and Lucius [Ref. 7] used numerical methods to obtain 
asymptotic solutions to the boundary layer equations for 
oscillating flow on a semi-infinite flat plate with no 
pressure gradient. Their solution agrees with those of 
Lighthill and Nickerson. 

The previously mentioned studies have all dealt with 
laminar boundary layer flows . The analytic approach to tur- 
bulent boundary layer flow has apparently been less fruitful. 
Karlssen [Ref. 8] separated the fluctuations into periodic 
and random components, and by averaging over a complete 
period obtained equations similar to the steady turbulent 
boundary layer equations, but with an additional fluctuating 
term. Other analytic works [Ref. 9, 10] have used linearized 
solutions to solve the pressure distribution on an airfoil 
at low angle of attack, but are applicable to only a small 
number of cases. 

2 . Experimental Studies 

A great deal more work has been done experimentally 
than analytically, but there are still large gaps in the 
field. Nickerson [Ref. 6] partially verified both his 
analysis and that of Lighthill. Miller [Ref. 11] studied • 
the transition phenomena on a flat plate in oscillating flow. 
He was able to confirm some of the previous analytic pre- 
dictions. In addition he determined the transition Reynolds 
number, turbulent Reynolds number, and turbulent intermit- 
tency factor for oscillating Blasius type flow. Despard 
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[Ref. 12] investigated the separation of a laminar boundary 
layer in oscillating flow. He proposed that the definition 
of separation for a laminar boundary layer in an oscillating 
freestream be the initial occurrence of zero velocity or 
reverse flow at some point in the velocity profile through- 
out the entire cycle of oscillation. He was also able to 
make some prediction about the behavior of the separation 
point . 

Morrissey [Ref. 13] investigated the effects of 
large amplitude flow oscillations on the heat transfer from 
a flat plate with a turbulent boundary layer. Jacobs [Ref. 
14] studied the effect of oscillating mean freestream on 
the turbulent intensity distribution in a turbulent boundary 
layer. These tvjo studies indicate that unsteady flows ex- 
hibit no significant alteration in the character of the eddy 
diffusivity or the turbulent intensity distribution when 
compared to a steady turbulent flow. Banning [Ref. 15] in- 
vestigated the pressure distribution on an airfoil in a 
turbulent oscillating freestream. Others have studied the 
effects of oscillating flows with compressibility effects 
included . 

B. BLOWING TO AVOID SEPARATION 
1 . Experimental Studies 

Boundary layer control through blowing has been with 
us almost as long as the boundary layer concept itself. The 
idea arose initially from the use of slots as a ].ift augmen- 
tation device, v.’hich was suggested by Lac}imann of Germany in 
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1918 and later tried by him and Handley-Page in England. 

Betz [Ref. 16] theorized that the effect of the slots was 
to accelerate the boundary layer. Baumann [Ref. 17] was 
led by this inter'pretation to replace the air passing through 
the slot with air ejected from the wing interior. This pro- 
duced roughly the same effect as the slots with the added 
advantage of being able to control the eA'ect by control of 
the pressure inside the wing. 

Until the 1940 's little more was done with blowing 
except some experiments which used it as an alternative for 
ailerons. This was done by blowing out a slot over a short 
wing section to induce high lift. By differential blowing 
on the wing a rolling moment was produced similar to that 
of an aileron. In 1942-43, triggered, no doubt, by the war, 
experiments began to be conducted in France, Germany and the 
United States, followed shortly by Britain. Probably the 
most important contribution of this era was made in 1948 by 
Poisson-Quinton [Ref. 18]. Prior to this time the common 
parameter to measure blowing was 



where is the volumetric flow through the slot, U^ the 
mean freestream velocity and S the characteristic area of 
interest. Schwier [Ref. 19] had shown that narrow slots 
were more efficient than v;ide ones. Poisson-Quinton verified 
this and suggested that a more appropriate blowing parameter 
was 




(II-l) 
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(II-2) 



m , V . 

C = - J J 

q«.s 

v;here m. is the mass flow of blown air, V. the iet velocity 
at the exit assuming isentropic expansion from the total 
pressure and temperature inside the jet supply duct, and 
the freestream dynamic pressure. Subsequent experimental 
v7ork supported the use of this parameter %o predict the 
effectiveness of blowing. 

Unfortunately all the early efforts in boundary 
layer control through blowing were hampei'ed by the great 
difficulties involved in practical applications. Major 
obstacles vjere a supply of high pressure air and the weight 
penalty associated with the ducting. However, with the 
advent of gas turbine engines for aircraft propulsion a 
good supply of air became available. With this development 
and advances in metal alloys and other lightweight materials 
of high strength, blowing has become much more attractive. 
Since 1950 the use of blowing has proceeded at an accelerat- 
ing pace. Unfortunately little work has been accomplished 
in the area with which this study is concerned. The majority 
of the v7ork has been in the use of blowing for lift augmenta- 
tion by delaying stall or by providing circulation control. 
All of these experiments present results in terms of 

and C and their variation with C . Also the values of C, 

m y VI • 

required for lift augmentation are much greater than that 
required to keep the boundary layer from separating under 
the influence of an adverse pressure gradient. As a result 
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there has been little interest in studying the amount of 
blowing required to maintain streamline flow. The only in- 
formation to be gleaned from this work was that higher pres- 



sure gradients require higher to overcome separation. 

2 . Analytic Studies 

Theoretically, most of the interest lies in increas- 
ing lift , rather than detailed study of blowing requirements 
to maintain attachment. There are, however, a few notable 
exceptions. Carriere and Eichelbrenner [Ref. 23] developed 
a method for calculating the conditions for flow re-attach- 
ment by a jet discharging against adverse pressure gradient. 
Their analysis, unfortunately, was heavily dependent on the 
availability of empirical profiles from experimentation and 
several arbitrary assumpti o^^s . Within these ].imits it did, 
however, provide solution to cases with weak pressure gra- 
dients. Kozlous and Zyanyak [Ref. 24] have recently anal- 
yzed a laminar boundary layer in an unsteady incoming motion 
around a body of arbitrary shape with either suction or 
blowing. They used a six degree polynomial to describe the 
boundary layer under these conditions . Their result was an 
integral equation which could be integrated directly, or 
reduced to a quadrature, depending on the velocity profile. 
They were able to verify their analysis for a symmetric 
wing in a start-up situation. This can be extended to an 
oscillating flow by varying the description of the incoming 
motion. Although derived for lift and drag predictions in 
laminar flows with blowing, this analysis shov/s great promise 
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for being extended to turbulent boundary layer blowing re- 
quirement predictions. The greatest difficulty with this 
extension, or any prediction of blowing in turbulent 
regimes , is that much of the energy imparted to the boun- 
dary layer is by turbulent mixing which defies accurate 
modeling . 

C. APPLICATIONS 

Recently there has been work accomplished toward the 
application of boundary layer control in an oscillating 
flow. Englar and Williams [Ref. 12, 23, 24] have applied 
boundary layer energization through blowing to augment the 
lift of a submarine stern plane and a symmetric airfoil at 
angle of attack. They have also applied tangential blovjing 
to the blades of a helicopter rotor system (the circulation 
control rotor) with encouraging results. 

There has not, however, been any detailed studies of 
the actual flow over the surface and its behavior under the 
conditions of high angle of attack, large pressure gradient, 
turbulence, oscillating freestream, or blowing to energize 
the boundary layer. The prediction of blowing requirements 
under these conditions will for the foreseeable future, 
rely heavily on empirical methods based upon experimental 
results. The present study is an attempt to provide some 
of this experimental data. 
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III. DESCRIPTION OF APPARATUS 



A. OSCILLATING FLOW WIND TUNNEL 
I . General Description 

The experimental work was conducted in the lowspeed, 
oscillating flow wind tunnel located in the Aeronautics 
Laboratories of the Naval Postgraduate School. This wind 
tunnel is of open circuit design, with a 24-inch square by 
223-inch long test section. A plan view of the tunnel is 
presented in Figure III.l. The tunnel inlet is eight feet 
square, resulting in a 16:1 contraction ratio. Three high 
solidity screens located in the inlet section just upstream 
of the nozzle produce measured freestream turbulence in- 
tensities of 0.261 to 0.413 percent for the velocities 
encountered in the present work. 

The wind tunnel drive consists of two Joy Axivane 
Fans in series, each of which has an internal, 100 horse- 
power, direct connected, 17.S0 rpm motor. The fan blades 
are internally adjustable through a pitch range of 25 to 
55 degrees, providing a v/ide operating range of test section 
velocities . Two sets of variable inlet vanes , located im- 
mediately upstream of each fan , are externally operated to 
provide control of test section velocity. These vanes are 
of multileaf design, and preswirl the air in the direction 
of fan rotation to reduce fan capacity. The total range of 
tunnel velocity is from 10 to 250 ft. per second. 
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2 . Rotating Shutter Valve 



Two fundamental methods of creating an oscillating 
flow environment have been employed in the past. Nickerson 
[Ref. 6] introduced oscillations by oscillating the model 
in a steady flow environment. This method severely re- 
stricts the range of attainable frequencies because of 
mechanical complications , and also introduces measurement 
difficulties . The other method is to actually oscillate the 
flow over a stationary model. Hill [Ref. 25] used a sliding 
shutter' to impose oscillations on the freestream but was re- 
stricted by mechanical limitations to low frequencies. 

The most successful method of obtaining an oscillat- 
ing flow with large ranges of frequency and amplitude was 
that employed by Karlssen [Ref. 8] and later by Miller [Ref. 
11] in his investigation of transition phenomena. A rotating 
shutter valve, immediately downstream of the test section, 
is employed to superimpose a periodic variation of velocity 
on the mean flow. The method used in the present investiga- 
tion is identical to that employed by Miller. The shutter 
valve consists of four horizontal steel shafts equally spaced 
across the test section. The shafts are slotted to accommo- 
date flat blades of various widths , forming a set of four • 
butterfly valves spanning the test section. Figure III. 2 is 
a photogx"*aph of the shutter valve . Each blade is driven from 
its immediate neighbor by means of a timing belt and pulley 
arrangement. The bottom shaft is driven by a five horse- 
power, variable speed, electric motor, througli a timing belt 



20 



and pulley. An intermediate shaft between the motor and 
shutter valve permits a wide variety of pulley ratios. This 
drive arrangement provides a frequency range from 0.4 cycles 
per second to the first critical frequency of 933 cycles 
per second. The electric motor presently in use, however, 
restricts the oscillation frequency to a maximum of 240 
cycles per second. The amplitude of oscillation is controlled 
by blade width. Test section closure may be varied from 25 
to 100 percent. The resulting amplitude of oscillation of 
test section velocity is a function of frequency, mean 
velocity, and pressure gradient. In this investigation 
blades producing 83.3 percent closure were used resulting in 
a perturbation range from 5 to 25 percent of the local mean 
freestream velocity. Figure III . 3 is a picture of the shut- 
ter valve drive arrangement. 

3 . Test Section 

The wind tunnel test section is shown in Figure III. 4. 
Continuous pieces of two-inch thick, aluminum, 24 inches wide 
and 223 inches long, form the upper and lower test section 
walls. Each of the sidevjalls consists of three two-inch 
thick panels of stress relieved lucite. For this investiga- 
tion the central sidevjall panel on the opposite side of the 
tunnel from the control console was replaced with two-inch 
thick plyvcood to facilitate the mounting of instrumentation 
and access to the model plenum. The Lucite panels on the 
console side of the test section are hinged and may be 
raised hydraulically, providing ac.cess to the test section. 
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Figure III. 5 shows the test section with the door open. 

The heavy construction of the test section is intended to 
minimize deflections induced by rapid changes in static 
pressure . 

Previous test on the test section velocity profile 
have shown the velocity variation is less than one percent 
of mean to within three inches of the wall [Ref. 12]. 

Figure III. 6 is a photographic viev/ of the test section and 
control console with the model in place. 

B. MODEL 

The model used in this investigation consisted of two 
plugs placed in the test section with a flat plate halfway 
between them. Each plug had a smoothly curved leading edge, 
followed by 42-inch straight section, then a diverging sec- 
tion as shown in Fig. III. 7. The maximum thickness of each 
plug was six inches and the plate was \ inch leaving 5-3/4 
inches between each plug and the plate. The diverging sec- 
tion departs from the horizontal by 27°. The slant length 
was 13 inches and the characteristic area chosen for calcu- 
lation of blowing coefficient was the plane area of the 
diverging section which was 305.5 square inches. 

Each plug contained a 1700 cubic inch plenum chamber. 

The entrance to the plenum was a two-inch diameter' hole to 
which the blovjing air supply hose was connected. The blow- 
ing slots were 2 3.5 inches vjide . The upper slot height v;as 
.041 inch and the lov?er v/as .055 inch. The difference in 
slot heights provided information on the blov/ing coefficient 
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requirements as a function of slot height. Because of its 
smaller slot height the upper jet would need less mass flow 
but higher jet velocity to achieve the same blowing coeffi- 
cient as the lower slot. Figure III. 8 shows the plenum and 
slot configuration. Blowing air was provided by a Carrier 
three-stage centrifugal compressor, driven by a 300 horse- 
power General Electric induction motor. It was capable of 
supplying 1900 cubic feet per minute at 29.5 psia. Each 
plenum was supplied and metered independently by a gate valve 
in the three-inch supply pipe. Each pipe had an orifice 
plate with a 1.8 inch hole diameter for use in mass flow 
measurements . 

The flat plate, constructed of a slab of phenolic material 
24-inches wide, 60-inches long, and \ inch thick had a rounded 
leading edge in order to avoid leading edge separation. The 
last five inches was tapered and hinged. For this study the 
hinged portion remained at zero angle of incidence with re- 
spect to the flat plate. The plate was mounted at the center 
line of the test section and ran from a point 7.5-inches from 
the leading edge of the plugs to a point five-inches back to 
the end of the plugs as shown in Fig. III. 7. 

Figure III. 9 shows a schematic of the test setup. 

C. INSTRUMENTATION AND CALIBRATION 

1 . Freestream Sensors 

A conventional pitot static tube located 22.5 inches 
aft of the model leading edge and midway between the flat 
plate and the upper model section was used to measure mean 
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freestream dynamic pressure. A hot wire probe 29.5 inches 
aft of the model leading edge and midway between the plate 
and the lower model section was used to determine magnitude 
of velocity perturbations during oscillations. The probe 
was connected to a locally .manufactured , single channel, 
hot wire anemometer which was connected to a Tektronix 555 
Dual-beam oscilloscope for display and read out. Previous 
experimentation [Ref. 12] has indicated that the hot wire 
circuitry is linear with freestream velocity. A short 
calibration run was made which verified this. 

The frequency of oscillation was obtained from a 
magnetic pickup mounted out-board of the uppermost shutter 
blade shaft, as seen in Fig. III. 2. The output frequency 
was read on a Dynascience digital counter. 

The pressure distribution in the tunnel was measured 
by a series of flush pressure ports one inch apart in the 
flat plate. They ran from a point five inches in front of 
the diverging section of the model to 2.3 inches from the 
end of the model. The ports were connected to a 40 tube, 

250 centimeter upright manometer board. The median reading 
of the oscillations was the recorded pressure value. 

2 . Blowing Coefficient 

The mass flow of blowing air vjas measured by a pair 
of orifice plates in the three inch supply lines. The ori- 
fice plates were calibrated in situ. This was accomplished 
by connecting the blov7er supply lines to a 76.5 cubic feet 
per minute rotameter manufactured by Fischer and Porter 



Company. The pressure drop across the plates was measured 
by a 36 inch U-tube water manometer. The results of the 
calibration were plotted against the ASME values as shown 
in Fig. III. 10. The differences were so small as to be 
negligible, therefore it was deemed reasonable to extrapolate 
the curve along the same form as the ASME curve to cover 
values not covered by the rotameter. Figure III. 10 was then 
used to find the volume flow of air for any value of pres- 
sure drop. 

The pressure downstream of the orifice plate was 
measured v;ith a 36 inch U-tube mercurcy manometer vented to 
the atmosphere. Air temperature remained essentially con- 
stant from the compressor to the plenums. With the volume 
flow, pressure and temperature, the mass flow could easily 
be calculated. 

The blowing jet velocity for calculation of has, 
historically, been defined as the velocity at the exit as- 
suming isentropic from some pressure and temperature inside 
the supply duct. These were measured with a special purpose 
Kiel temperature probe. The pressure output was connected 
to the same upright manometer board as the pressure distri- 
bution ports. The thermocouple for temperature measurement 
was connected to a Leeds and Northrup portable potentiometer. 

3 . Visualization of Separation 

Separation lines v/ere visualized with soft twine 
tufts placed along the diverging planes. The lines of tufts 
spanned the test section to within three inclics of each 
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sidev;all and were spaced at one inch intervals across the 
test section and two inch intervals down the plane as seen 
in Fig. III. 11. 
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Figure III.l. Plan View of Wind Tunnel. 
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Photograph of Rotating Shutter Valve. 





Figure III. 3. Photograph of Shutter Valve Drive Mechanism. 
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Figure III. 4. Photograph of Wind Tunnel Test Section Closed 
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Figure III. 5. Photograph of Test Section Open with Model. 
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Figure III. 6. Photograph of Test Section and Control Console. 
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Figure III. 7. Drawing of Model and Plate. 
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Figure III. 8. Drav;ing of Plenum and Blowing Slot Configuration, 



Plenum Temperature 
and Pressure Probe 
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Figure III. 9. Schematic of Test Setup. 
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Figure III. 10. Pressure Drop Across Orifice vs. Volume Flow of Air. 







Figure III. 11. Photograph of Diverging Planes from Downstream. 
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IV. TEST PROGRAM 



A. EXPERIMENTAL PROCEDURE 

The basic test program was to find, for some value of 
freestream dynamic pressure, the value of blowing coefficient 
that was just sufficient to maintain an attached boundary 
layer, then to study the effects of varying the frequency 
of oscillation on the blowing coefficient required. 

The first problem v/as to select values of mean free- 
stream dynamic pressure for the study. The dynamic pres- 
sures chosen were those which indicated 5, 7.5, and 10 cen- 
timeters on the micro-manometer. This decision was based 
on three criteria: 

1. Below 5 centimeters the atmospheric wind conditions 
at the inlet began to affect the test section conditions. 

2. Above 10 centimeters the fans could not maintain the 
pressure at higher oscillation frequencies. This could be 
remedied by opening the fan housing and altering the vane 
direction but this would have vastly increased the time re- 
quired for runs and would have had an effect on test section 
conditions . 

3. These three settings were easily read and maintained 
during tunnel operations. 

These three dynamic pressures gave values of mean fi'eestream 
velocity of about 96, 118, and 138 feet per second, respec- 
tively, although this varied slightly witli oscillation 
frequency . 
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The initial tests were dictated by the wind tunnel set- 
up from previous work. The set of belts and pulleys con- 
necting the electric drive motor to the shutter valve, at 
the beginning, corresponded to a frequency range of 2.48 to 
14.4 cycles per second. Then on the basis of the previous 
test results the range was expanded outward to cover fre- 
quencies of 1.6 to 16 cycles per second. (See Results for 
detailed explanation.) 

The detailed test procedure was as follows: 

1. Start tunnel and stabilize at desired dynamic 
pressure . 

2. Adjust blowing to barely maintain attached flow 
down the diverging plane. 

3. Take data for pressure distribution, blower plenum 
conditions and mass flow of blown air. 

4. Start shutter valve and set to the desired frequency. 
Adjust dynamic pressure as necessary. 

5 . Repeat step 2 . 

6. Take data for magnitude of oscillation. 

7. Repeat step 3. 

8. Adjust shutter valve for new frequency and repeat 
steps 5, 6, and 7 in order. 

B. DATA REDUCTION 

The data reduction v/as accomplished using a Fortran 
language computer program on the IBM 360/67 computer located 
at the .Naval Postgraduate School. The program is shown in 
a separate section beginning on page 61. 
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The calculations involved are of mean freestream static 
pressure, density, velocity and magnitude of the oscillation 
of mean freestream velocity in percent. Then the pressure 
distribution was non-dimensionalized with respect to free- 
stream dynamic pressure. Finally the blowing coefficient 
was calculated for the upper and lower blowers. 

In addition the computer program plotted the pressure 
distribution for constant frequency and blowing coefficient 
versus frequency for constant freestream dynamic pressure. 
The pressure distribution plots were used for further, manua 
data reduction. The blowing coefficient variation with fre- 
quency was a guide to further experimentation in the early 
stages, and for presentation of results. 



V. RESULTS 



The results of this study are presented in Figures V.2 
through V.ll. 

One of the greatest difficulties encountered during ex- 
perimentation was establishing a reasonable definition of 
separation in an oscillating flow. Turbulent oscillating 
flow does not lend itself to a single, unambiguous defini- 
tion of separation because the flow tends to separate pro- 
gressively. Flow areas may be identified that are attached 
always , separated always , and alternately attached and 
separated during a single cycle of oscillation. The method 
of visualizing separation with rows of tufts led to the 
definition used in this study. It was assumed that separa- 
tion had occurred if the flow reversed itself within the 
boundary layer during any part of a cycle. The definition 
used in this study does not coincide with Despard's defini- 
tion [Ref. 12] for separation in laminar boundary layers. 
Despard proposed separation as commencing with the initial 
occurrence of zero velocity or reverse flow throughout the 
entire cycle. Despard's definition cannot be easily used 
for flov7 visualization studies. Preliminary testing re- 
vealed that below about 1.6 cycles per second the flov; acted 
in a quasi-steady manner. That is, the blowing required to 
maintain attachment was dependent on the instantaneous 
velocity. Therefore it was necessary to provide blowing 
equal to that required by the maximum magnitude of the advei’se 
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pressure gradient during a cycle. Also, above 16 cycles 
per second, the oscillations were so rapid that it became 
impossible to detect, by eye, the onset of separation with 
the method in use. In fact, as the frequency approached 16 
cycles per second the decision as to whether to flow was or 
was not separated became more arbitrary because the tufts 
would oscillate at a frequency close to t^*ie limit of the 
eyes ability to discern change. For the reasons cited above 
the study was continued to the frequency range from 1.6 to 
16 cycles per second. 

One other problem encountered with the experimental setup 
was that of three dimensional flow in the test section. It 
has been suggested that the study of separation on a flat 
plate in a turbulent oscillating flow with sLi'ong advei'se 
pressure gradient could be accomplished by using a model 
similar to that used in this study. To achieve this it 
would be necessary to have the flow attached across the 
entire plane of the diverging section. This should allow 
the boundary layer growing on the flat plate to become de- 
tached due to the impressed pressure gradient. Tufts were 
placed, therefore, on both the upper and lower surfaces of 
the flat plate, as well as along the diverging plane. Boun- 
dary layer growth in the corners and on the sidewalls caused 
a wedge shaped area of separation to appear on the divei^ging 
plane as shown in Figure V.l. Tufts close to the sidewalls 
indicated this corner-sidewall separation but those close to 
the centerline of the tunnel indicated the flow v;as still 



attached to the diverging plane. It was necessary, therefore, 
to observe only tufts in the center of the test section dur- 
ing this study. The flow was considered attached if at least 
five contiguous tufts in one row indicated that the flow was 
attached. This definition was somewhat arbitrary, but its 
consistent application led to some meaningful trends . Some 
separation did occur from the bottom side of the plate when 
the mean freestream velocity was 138 feet per second with 
high blowing coefficient. Unfortunately, this separation 
was not reproducible. The flow never did detach from the 
top side of the plate. These observations show that it will 
be necessary to blow along the sidewalls to prevent the cor- 
ner separation, thereby ensuring more predictable separation 
characteristics on the plate. 

To analyze the results it was necessary to determine the 
parameters on which C (the minimum blowing required to 
maintain attached flov;) is dependent. It was assumed that 
the jet velocity required V. to attach the flow was a func- 
tion of the density, p; mean freestream flow velocity, ; 
the length of the plane on which attachment is measured, £; 
viscosity, y; the magnitude of the perturbation of the free- 
stream velocity, AU^ ; the frequency of the oscillation, co ; 
the slot height, h; and the effect of the pressure gradient 
as measured by APj,*6, v/here 6 is the boundary layer thickness 
at the diffuser exit and AP^ the pressure difference from 
the reference or undisturbed condition to the exit. This 
dependence is expressed as Equation V-1. 
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(V-1) 



V. = f.[p,U^,£,AP •6,y,AU^,w,h]. 

Jr J- L 

Bernoulli’s equation may be used to express the pressure 
coefficient as 



El = r 1 _ _E 1 

?5PU2 U2 



(V-2) 



where is the velocity at the exit of the diverging section. 
Substituting V-2 into V-1 and non-dimensionalizing with re- 
spect to p, U^, and £ leads to 

(V-3) 
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It is convenient to replace y/pU^£ by its reciprocal, the 
Reynolds number. 

C may then be expressed as 



y 
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ih.V. pbhV. 
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where b is the width of the slot. This may be simplified to 

V. 



C = 2 ( — H ^ 



(V-5) 



From equation V-3 and V-5 the critical blowing coefficient 
may be v.’ritten as 

AU 
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Equation V-6 shov/s that C is a function of five independent 

’■'r 

parameters. In the experimental study the pressure gradient 



parameter 



. Up 

r ^ ) r 1 _ _±i I 

00 

was determined by the model and tunnel configuration. Al- 
though the pressure gradient was measured (see later dis- 
cussion) there was no measurement or control of 6 possible 
in this study. This parameter may be most important and 

further studies should be made while observing 6 . Three 

5 5 

Reynolds numbers were considered; 5.93x10 , 7.31x10 and 

5 

8.53x10 . These Reynolds numbers corresponded to test sec- 
tion freestream velocities of 96 fps , 118 fps, and 138 feet 
per second (dynamic pressures of 5.0, 7.5, and 10.0 centi- 
meters of water), respectively. 

Oscillation frequency was varied continuously between 
1.6 and 16 cycles per second which resulted in a reduced 
frequency range of approximately .1 to 1.1. The percent 
perturbation of velocity AU^/U^ was a function of frequency, 
shutter valve plate size, and freestream dynamic pressure. 

It was not independently controllable in this study since a 
single plate size was used. therefore was dependent 

on frequency and dynamic pressure. Since AU^/U^ did vary as 
the reduced frequency was changed, the resulting frequency 

variation in C measured could not be determined at constant 
^R 

values of all the correlation parameters. 

Reference [18] states that experiments have shown that 
for steady flows C is not a function of h/Z for high slot 
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Reynolds Numbers, V^h/v. Therefore, Y\/Z is normally omitted 
from Equation V-6 . Blov.’ing characteristics could vary sig- 
nificantly with slot height at lov? values of the slot Reynolds 

4 5 



numbers. Slot Reynolds numbers encountered in this study varied 

from 6.06x10^ to 11.56x10^ for the upper slot, and 6,42x10^ to 
3 

11.74x10 for the lower slot. To independently check the re- 
sults of Ref. [18] for the case of unsteady flows, this study 

- 3 - 3 

used two non-dimensional slot heights of 3.15x10 and 4.23x10 
Figures V-2, 3 and 4 clearly indicate that frequency has 
a definite effect on blowing requirements to maintain an at- 
tached flow. At a Reynolds number of 5.93x10^, Figure V-2 
shows that there is no apparent change in C from the non- 
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oscillating condition to a reduced frequency of 0.5 where 

C increases by approximately 30 percent. The required 
^R 

flowing gradually decreases through the remaindei’ of the 
frequencies tested in this study. For the Reynolds number 

5 

of 7.3]xl0 . Figure V-3, the initial values of C at a 

4 u 

reduced frequency of 0.1 are approximately 40 percent 

above the non-oscillating condition. At a reduced frequency 

of 0.13, C drops to approximately 30 percent greater. At 
^R 

0.48 the blowing required increases to as much as 50 percent 
above the non-oscillating condition after which it drops 
back down to approximately 30 percent for the higher fre- 
quency oscillations. Figure V-4 shows the C behavior for 
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a Reynolds number of 8.53x10 . It shows that a reduced 

frequency of 0.08, C is approximately 40 percent greater 

^R 

than the non-oscillating case. C rises rapidly to 100 

percent higher at a reduced frequency of 0.12. C stays at 

^R 

this level until a reduced frequency of 0.2 9 is reaclied 
vr/here it drops to 50 percent above the ]ion-oscillating 
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condition. At a reduced frequency of 0.4, the blowing re- 
quired then rises rapidly to 70 percent above the non-oscil- 
lating condition followed by a gradual decrease through the 
remainder of the frequencies studied. 

There was a significant difference between C for the 
upper and lower blowers. This was not expected, in spite 
of the differences in slot heights, as indicated in pre- 
vious discussion, because is a momentum coefficient. 
Therefore, although the upper blower had less mass flow, 
its greater jet velocity should produce the same attachment 
effectiveness as the lower blower. The upper blov;er C 

^ R 

was approximately 15 percent higher than that of the lower 
blower at R^^ of 5.93x10^. For R^ of 8.53x10^ the difference 
was 20 percenr. But at R^ equal to 7.31x10^, the upper 
blovjer required approximately 4 percent less blowing than 
the lower blower. This inconsistency, with test section 
velocity, of the difference between upper and lower C 

^R 

might indicate that the upper and lower diverging sections 

experience different flow environments. Since separation 

is intimately linked with both previous boundary layer 

growth and disturbing influences, the C differences may 

^R 

be explained by the presence of some downstream flow asym- 
metry that is a function of test section velocity. 

It was recognized that the variation in C was depen- 

R 

dent on the pressure gradient and the magnitude of the 
perturbation of freestream velocity, both of which exhibit 
variations vjith frequency in this test £’.etup. The median 
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pressure coefficient, , on the flat plate was plotted versus 

position, X, along the plate. x/i, equal to zero corresponds 

to the blowing slot location. Typical pressure distributions 

are shown in Figure V-5 for the cases of steady flow with no 

blowing, steady flow with blowing, and oscillating flow with 

blowing respectively. The pressure gradient was determined 

by estimating the slope of the best straight line fit of the 

plotted pressures. The resulting slope was plotted against 

the reduced frequency for each of the freestream velocities 

run. These data showed similar variation of C versus re- 

^x 

duced frequency and so are plotted together in Figure V-6. 

An examtination of Figure V-6 shows that the pressure gradient 
variation with frequency was not the primary cause of the 
exhibited large variations of C with fi'equenuy. Similar 
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plots of percent perturbation in freestream velocity, as 
measured by the hot wire, versus reduced frequency are shown 

5 

in Figures V-7, 8 and 9. The cases of equal to 5.93x10 
and 7.31x10^ show a gradual, almost linear, decrease of 
perturbation velocity with frequency but the R^^^ of 8.5 3x10^ 
case shovjed the characteristics seen in Figure V-9. There 
was, however, no indication that these variations in the 
velocity perturbation were responsible for the shape of the 
large variation of C with frequency seen in Figure V-2 , 3 

’■‘r 

and 4 . 

The blov;ing requirement to maintain attached flow has 
been shown to be frequency dependent. This dependence ex- 
hibits characteristics whicli sue’gost resonant behavior. 
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This resonance might be the result of flow confinement. 

Wave reflections, due to the tunrel walls, could affect the 
conditions in the test section. These reflected waves could 
either reinforce or dissipate the shutter valve induced 
pressure pulses. Another source of resonance lies in the 
physical makeup of the tunnel. It was observed during the 
experimental runs that certain frequency bands at each tun- 
nel test section velocity made the wind tunnel vibrate a 
great deal, sometimes causing the test section to be travel- 
ing as much as two inches longitudinally. These resonant 
frequency bands are shown in Figures V-2, 3 and 4. This 
resonant behavior of the tunnel occurred around a frequency 
of approximately 8.7 cps. Tunnel vibrational modes are seen 
to be influencing the blowing characteristics observed. Dif- 
ferent test configurations may exhibit different characteristics. 
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Figure V.l. Photograph of Tuft Showing Corner Sidewall 
Separation. 
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Figure V.2. Blowing Coefficient Required vs. Frequency, 
Rjj^ = 5.9 3x10^. 
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Figure V.3. 



Blov;ing Coefficient Required vs. Frequency, 
R^ = 7.31x10^. 
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Figure V.4. Blov.’ing Coefficient Required vs. Frequency, 
R^ = 8.53x105. 
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Figure V.5. Pressure Coefficient vs. Position. 
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Figure V.6. Pressure Gradient vs. Frequency 
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.7. Velocity Perturbation vs. Frequency, R = 5.93x10^. 
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8. Velocity Perturbation vs. Frequency, R^ = 7.31x10' 
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Figure V.9. Velocity Perturbation vs. Frequency > 
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VI. CONCLUSIONS 



From the results obtained, the following conclusions may 
be drawn: 

1. There is a definite effect on the blowing required 
to maintain an attached boundary layer in a strong adverse 
pressure gradient as oscillations are superimposed on the 
freestream velocity and the frequency is varied. 

2 . The frequency dependence of the blowing requirements 
exhibits characteristics which suggest resonant behavior, 

3. Considerable three-dimensional flow is produced in 
the present test setup and sidewall blowing will be required 
to produce two-dimensional separation on the plate. 
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VII. REC OMMEWI*/' ON S 



1. C is a function of the five parameters discussed 



y 



R 



in the Results section. More tests should be made in which 
only one of these is varied at a time. As discussed earlier, 
measurement of the boundary layer thickness is most desirable. 

2 . The velocity profiles throughout the flow field need 
to be studied as a function of time while varying oscillation 
frequency. To accomplish this a method for flow measurement 
must be developed which will fulfill the following require- 
ments; non-interference with the flow, capable of measuring 
turbulent variations , capable of scanning a relatively large 
section of the flow, and capable of resolving rapid changes 
in flow velocities at high oscillation frequencies. 

3. An investigation of the natural frequencies of the 
tunnel-model-air flow configuration should be made to deter- 
mine the effects on pressure gradient, magnitude of oscil- 
lating velocity perturbation, and blowing requirements to 
maintain an attached boundary layer. 

4. Sidewall blowing should be used to ensure two dimen- 
sional separation of the flat plate for further studies. 
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COMPUTER PROGRAM 



Q :{:3!c :!!<c *s!t :{:5!: « >!i sJc :{c :Jc ait jjt ^5}: sjt jjtsft iji sisJ^! ;|t s{t s{t 

C « 

C « 

C * THIS PROGRAM CALCULATES THE FOLLOWING . 

'' "• 1. MEAN FREESTREAM VELOCITY IN FEET PER SECOND 

2. OSCILLATION IN RADIANS PER SECOND 

3. FREESTREAM VELOCITY PERTLRBITION IN PERCENT 

4. DIMENSIONLESS PRESSURE DISTRIBUTION 
A. ALSO PLOTTED ON CALCCMP PLOTTER WITH THE 

ABSCISSA THE X POSITION IN INCHES 
BEGINNING AT THE UPSTREAM POSITION AND 
THE ORDINATE IS THE PRESSURE COEFFICIENT 

5. JET VELOCITY OF EJECTED AIR FOR BOTH UPPER 
AND LOWER BLOWERS IN FEET PER SECOND 

6. MASS FLOW OF EJECTED AIR IN 
POUNDS MASS/SECOND 

7. BLOWING COEFFICIENT FOR UPPER AND 
LOWER BLOWERS 

A. PLOTTED ON CALCOMP PLOTTER WITH THE 
ABSCISSA THE FREQUENCY IN RADIANS AND 
THE ORDINATE THE BLOWING COEFFICIENT 
OUTPUT IS LINE PRINTED, SEPARATED BY FREQUENCY 
AND IS LABELLED 
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INPUT REQUIRED 

PATM = ATMOSPHERIC PRESSURE (INCHES OF MERCURY) 

T = OUTSIDE AIR TEMPERATURE (DEGREES FAREHEIT) 

FREQ = FREQUENCY (CYCLES PER SECOND) 

Q = MEAN FREESTREAM DYNAMIC PRESSURE (CENTIMETERS 
OF WATER) 

QMAX = MAXIMUM MAGNITUDE OF VELOCITY OSCiLATIONS 
(VOLTS OUTPUT OF HOTWIRE ANEMOMETER) 

QMIN = MINIMUM MAGNITUDE OF VELOCITY OSCILATIONS 
(AS ABOVE) 

PAMB = AMBIENT PRESSURE (CENTIMETERS OF WATER, 

REFERENCE HEIGHT FROM MULTI-TUBE MANOMETER BOARD 
P = PRESSURE DISTRIBUTION, 16 VALUES PUT IN 1X16 
VECTOR (CENTIMETERS OF WATER) 

P2U = PRESSURE JUST DOWNSTREAM OF THE ORIFICE PLATE 
- IN THE SUPPLY LINE TO THE UPPER BLOWER 
(CENTIMETERS OF WATER) 

P2L = PRESSURE AS ABOVE BUT LOWER BLOWER 
VFLOU = VOLUME FLOW OF BLOWING AIR TO UPPER BLOWER 
(CUBIC FEET PER MINUTE) 

VFLOL = VOLUME FLOW AS ABOVE FOR LOWER BLOWER 
TPU = PLENUM AIR TEMPERATURE, UPPER BLOWER PLENUM 
(DEGREES FARENHEIT) 

TPL = PLENUM AIK TEMPERATURE, LOWER BLOWER 
PPU = PLENUM AIR PRESSURE, UPPER BLOWER. 

(CENTIMEPS OF WATER DIFFERENCE FROM ATMOSPHERIC) 
PPL = PLENUM AIR PRESSURE, LOWER BLOWER 

INPUT GOES ON FOUR DATA CARDS, THE FIRST HAS 7 NUMBERS 
THEY ARE, PATM, T, FRFQ, Q, QMAX, QMIN, FAME. 

THE NEXT TWO CARD HAVE THE 16 MUMBEPS OF PRESSURE 
D 1STRI3UT luN. THE LAST CARD HAS, P2U, VFLOL, TPU, 

PPU, P2L, VFLOL, TPL, PPL. THE NUMbtRS ARE INPUT 
AS F10.4 
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LINE OUTPUT 

FREQUENCY (CYCLES PER SECOND) 

Q (CENTIMETERS OF WATER) 

VELOCITY (MEAN FREESTREAM IN FEET PER SECOND) 
PERCENT (PERTURBATION OF VELOCITY) 

OMEGA (OSCILATICNS IN RADIANS PER SECOND) 

CMU (BLOWING COEFFICIENT, UPPER BLOWER) 

VJU (JET VELOCITY, UPPER BLOWER IN FEET/SECOND 
FLOU (MASS FLOW OF JET, UPPER BLOWER IN 
POUNDS PEP SECOND) 

CML (BLOWING COEFFICIENT, LOWER) 

VJL (JET VELOCITY, LOWER) 

FLOL (JET MASS FLOW, LOWER) 

STATIONS (PRESSURE D I STR I BUT I ON ^TAP POSITIONS 
WITH 1 AS THE DOWNSTREAM POSITION) 
CP (PRESSURE COEFFICENT AT STATION) 



PLOTTED OUTPUT 

PRESSURE COEFFICIENT VS X POSITION (WITH 1 BEING 
UPSTREAM) 

BLOWING COEFFICIENT VS OMEGA (PLOT OF CML AND CML 
VERSUS OMEGA FOR EACH INPUT DECK) 



DIMENSION CP(16,40),CMU(40),OMEGA(40),F(40) , 
1PCSL(16),P(16),CML(40) ,LPOS(16) , C P 1 ( 1 6 ) , QA ( 40 ) 

10 FORMAT (7F10.4) 

15 FORMAT (3F10.4) 

20 FORMAT ( • 0 ' , T3 , * Q • , T 1 7 , • VELOC I T Y » , T3 0 , » P E RC ENT ' , T44 , 
I'OMEGA* ,T58, 'CMU' , T72 , 'VJU« ,T84, 'FLOU' ,T98, 'CML' 
i,T110, 'VJL' ,T123, *FLOL‘ ) 

25 FORMAT ( • 0 ' , lOF 13 . 5) 

30 FORMAT ( ' 0 • , T 1 9 , • STAT I ON ' , T37 , ' C P • / ) 

35 FORMAT (' • ,I20,F26.8) 

40 FORMAT Cl') 

45 FORMAT (6A8) 

50 FORMAT COLAST = »,I10) 

55 FORMAT COTHE FREQUENCY IS ',F10.5) 

REAL--:^8 LABELl/'CML ’/ 

REALMS LA6EL2/'CMU •/ 

REAL>^8 LABELS/' •/ 

REALMS TITLEK12) 

READ (5,45) TITLEl 
REAL>!=8 TITLE2(12) 

READ (5,45) TITLE2 
DO 1000 J=l,28 

READ (5,10) PATM,T,FREQ,Q,OMAX,QMIN, PAMB 

READ (5,15) P 

DC 100 K=l,16 

LPOS(K) =K 

PCSL (K) =LPOS(K) 

100 CONTINUE 

READ (5,15) P2U,VFL0U,TPU,PPU,P2L,VFL0L,TPL,PPL 

CONVERT TEMPERATURES TO RANKINE 

TEMP=T+459.688 
TPL=TPL+459.688 
TPU=TPU+459.688 
QA( J) =Q 
F( J)=FREQ 

CALCULATE OMEGA, MEAN FREESTREAM VELOCITY, AND 
PERCENT PERTURBATION 

OMEGA (J )=FRE0*6.283185308 
PINF=^(P(14) + P(15)+P(16) )/3.0 
• RHO=( (2.0*PAMB) -PI NF) >'-PATM«.o41206962D0/PAME/TEMP 
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VEL=2.0 20916 9^4D0=«'SCRT (Q/RHO ) 

PERCEN= ( (QMAX-QMIN)/( GMAX+QMIM) )*100 

NCN DIMENSICNALIZe PRESSURE 

DC 200 1=1,16 
CP(I, J)=( PINF-P( I ) )/Q 
200 CONTINUE 

CALCULATE BLOWING PARAMETERS 

RH02U= (P2U<-PATM)*. 0412 06962 IDO/ TPU 
FLOU=VFLOU*RH02U/60.0 

VJU = DSQRT (5454. 5 0724D0*TPU^ (1.0-( ( ( 2. 0*PAMB I-PINF ) 
1/ ( (2.0'l=PAMB)-PPU ) ) *^.286224125600) ) 

CMU(J) = (FLGU*VJU)/ I 0*4.332271417 DOJ 
RH02L= (P2L+PATM)*.0412069621D0/TPL 
FLGL=VFL0L*RH02L/60.0 

VJL = 0SQRT(5454.5 0724D0*TPL* (i.0-( ( ( 2 . 0* P AMB ) - PI NF ) 
l/( (2.0*PAMB)-PPL ) )**. 286224 1256D0) ) 

CML(J)= (FL0L*VJL )/(Q*4. 3322714 17D0) 

LIKE OUTPUT 

WRITE (6,40) 

WRITE (6,55) F(J) 

WRITE (6,20) 

WRITE (6,25) Q,VEL, PERCEN,OMEGA( J),CMU( J) , 
1VJU,FL0U,Ci;l( J) ,VJL,FL0L 
WRITE (6,30) 

00 300 1=1,16 

WRITE (6,35) LPOS ( I ) ,CP(I , J) 

300 CONTINUE 
1000 CONTINUE 

PLOTTING ROUTINES 

NPTS=23 
EXSC=15.0 
YSCL=.006 
MX=0 
MY = 0 

CALL ORAW(NPTS,OMEGA,CML, 1, 1,LABEL1,TITLE1 
1 , EX SC , YSC L , 0 , 0 , MX , MY , 9 , 15 , 0 , LAS T ) 

WRITE (6,50) LAST 

CALL 0RAW(NPTS,DMEGA,CMU,3,2,LABEL2,TITLE1 
1,EXSC,YSCL,0,0,MX,MY,9,15,0,LAST) 

WRITE (6,50) LAST 
DO 2000 J=l,28 
DO 1100 1=1,16 
K=17-I 

CPI ( I ) =CP{ K, J) 

1100 CONTINUE 

CALL ORAW( 16, POSL ,CP1,0,2,LABEL3,TITLE2, 
12.0,0.0,0,0,i,0,8,8,0,LASl ) 

WRITE (6,50) LAST 
2000 CONTINUE 
STOP 
END 
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